Nonanalytic quantum oscillator image of complete replica symmetry breaking 



R. Oppermann, H. Schenck 

We describe the effect of replica symmetry breaking in the field distribution function P(h) of the 
T = SK-model as the difference between a split Gaussian and the first excited state ip\ of a weakly 
anharmonic oscillator with nonanalytic shift by means of the analogy P(h) o \tp 1 (x)\. 

New numerical calculations of the leading 100 orders of replica symmetry breaking (RSB) were 
performed in order to obtain P(h), employing the exact mapping between density of states p(E) of 
the fermionic SK-model and P(h) of the standard model, as derived by Perez-Castillo and Sherring- 
ton. Fast convergence towards a fixed point function p(E) for infinite steps of RSB is observed. A 
surprisingly small number of harmonic oscillator wave- functions suffices to represent this fixed point 
function. This allows to determine an anharmonic potential V(x) with nonanalytic shift, whose first 
excited state represents p(E) and hence P(h). The harmonic potential with unconventional shift 
V2(x) ~ (\x\ — xo) 2 = (x — xo sign(x)) 2 yields already a very good approximation, since anharmonic 
couplings of V(x) — V2(x) ~ \x\ m , m > 2, decay rapidly with increasing m. 
f*"*) ' We compare the pseudogap-forming effect of replica symmetry breaking, hosted by the fermionic 

SK-model, with the analogous effect in the Coulomb glass as designed by Davies-Lee-Rice and 
described by Miiller-Pankov. 
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I. INTRODUCTION 



The Sherrington-Kirkpatrick (SK-) models has been an extraordinary stimulating source for analytical and numer- 
ical studies in the physics of frustrated magnetism and in numerous interdisciplinary applications. Many relevant 
L '. properties of the SK-model have been explored since its discovery^—. 

Near the zero temperature limit and at T = 0, the model displays rather simple features but also analytically still 
unresolved behavior, and therefore remains challenging in our point of view. 

Along a particular line of theoretical analysis, on which the present article evolves, we searched for explicit analytical 
solutions deep in the ordered phase and particularly at zero temperature. The limit of zero temperature was recently 
described to accommodate two critical points related to the criticality of the hierarchical scheme of replica symmetry 
breaking (RSB) hosted by the SK-model^— . In order to reveal this type of criticality, the finite-step discrete-RSB 
was utilized. Extreme high order calculations up to 200 orders and Pade-approximants allowed to analyze the limit 
towards infinite breaking analyzed by means of high-precision numerical work, renormalization group ideas, and 
O ' scaling theory in addition^ - — . The scaling theory allowed to resolve the puzzle of noncommuting zero temperature 
limit and infinite-step limit of replica symmetry breaking^. 
' In this article we present new numerical calculations which led us to a surprisingly simple modeling of the distribution 
, of internal magnetic fields. 

Pioneering work on this local magnetic field distribution by Thomsen et al^ presented various analytical limits and 
numerical studies (TAP theory, simulations etcetera) for vector spin glasses (inch SK of course) and for the relevant 
temperature range from zero to above the glass transition. 

By virtue of an exact mapping as derived by Perez-Castillo and Sherrington^, we can, for our convenience, evaluate 
the spectral function (or density of states) for the fermionic version of the SK-model (see below for definition) and 
directly employ its functional form in one-to-one correspondence as the field distribution of the standard SK-model. 

In the following section we recall a few necessary details of the fermionic SK-model (sometimes called "ISG/" or 
" fSK" ) and of the fermionized SK-model as an exact representation of the standard SK-model. The formalism on 
how to deal with the Parisi scheme in the Grassmann field theories of these models has been presented in numerous 

■ earlier papers at length. 

Our main goal here is to report an astonishingly simple representation of the internal magnetic field distribution 

■ by the functional form of a single eigenstate of a weakly anharmonic quantum oscillator. 

This idea was guided by the known simple split- Gaussian replica-symmetric solution for the density of states p(E) 
at T = 0. The broad energy gap of this lowest approximation is filled by replica symmetry breaking. The resulting 
pseudogap reflects the criticality of RSB at T = 0. Its power law must be compared with the Efros-Shlovskii pseudogap, 
and indeed we will add below remarks on the Coulomb glass pseudogap as derived by Miiller and Pankov 1 -. Their 
intricate selfconsistent scheme 1 ^ is used to unravel the special role of the Coulomb interaction and of RSB for the 
pseudogap-creation at very low temperatures. 

The representation of the spectral density by means of ID quantum oscillator wavefunctions was also motivated by 
the advantage that one can integrate analytically the Lehmann representation of the single fermion Greens function 
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(term by term) for the fermionic SK-model. 
The paper is organized as follows: 

The main results on the functional mapping between field-distribution P(h), density of states p(E) and the first 
excited state of an unconventional quantum oscillator with nonanalytic shift and very weak anharrnonicity, is addressed 
in all chapters. The numerical data, on which the partially-analytical proposal is based, are presented and analyzed 
in section 4. The analytical formulas, which were evaluated numerically, are presented in section 11111 The theoretical 
basis which led to these analytical equations was published earlier— and employs the Grassmann-field and generating 
functionals so that the density of states can be extracted from the single-fermion propagator. Section 5 finally contains 
the comparison with the Coulomb glass. 



II. THE MODELS UNDER CONSIDERATION, AND THEIR RELATIONSHIPS 

A. The spin glass models 

Let us begin with the SK-model, its exact fermionic representation (fermionized SK-model), and the generic 
fermionic SK-model, sometimes called ISG/ of /SK-model. These models agree for half-filling (/i = 0) and at T = 0, 
and hence can be chosen for technical convenience. 

We took advantage of the grand-canonical apparatus of generating functionals in Grassmann field theory, which has 
been elaborated for the fermionized version. The fermionic SK-model, as the extreme localized limit of an itinerant 
spin glass model, contains certain aspects which can be compared with the Coulomb glass. 

The SK-model together with its grand-canonical fermionic representation (fcrmionization) and its fermionic variant, 
listed by the spin Hamiltonian Hsk and the fermionic partners T-LfSK respectively, are given by 

N 

Hsk = -'^JtjSiSj-'^HiSi with Si = ±l (1) 

i,j i=l 
N 

UfSK = - ^ ■l,,rr,T l - y^XHi °i - V n i) witn °» = n it ~ n ii ( 2 ) 

i,j i=l 

where Oi = — rijj., rii\ = a\ x ai\ in terms of the fermionic operators obeying {a a , a£} — S arr The spin interaction 
Jij is assumed to be range-free and Gaussian-distributed with zero mean (for simplicity) and variance J 2 /N. 

The chemical potential [a, assumed as nonrandom and homogeneous in this case, is chosen for the grand-canonical 
model Hamiltonian (2) in two ways as given by 

(2a) /i = i-^T fermionized SK-model 
(26) fi = half-filled SK-model 

The ingenious trick introduced by Popov and Fedotov^., which allows fermionization for /i = ^ttT, has been used 
before in many of our preceding papers (including its generalization to higher spin quantum numbers). This mapping 
is exact in the sense that the thermodynamics of the original model and of its fermionized version are identical. 

The 2nd model is also classical, since all operators commute - quantum time-dependence can only be observed when 
correlations of odd numbers of fermionic operators at different times are considered. The chemical potential signals 
that a grand-canonical ensemble is employed. For half-filling at T = 0, realized by /i = 0, both models agree. 

The distribution of internal magnetic fields P(h)& is defined for the SK-model by 

i 3 

while the fermionic density of states of the fermionic SK-model is defined via the disorder-averaged single fermion 
Greens function by 

p{E) = - E a )) = ( G*(E) ) (4) 

a 

where G R represents for the retarded and real-space-local fermionic Greens function. The bracket refers to statistical- 
and disorder-average over the Gaussian-distributed Jij . Standard many body formalism tells us how to obtain G R 
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from the real-time-ordered T — Greens function Gij(t — t') — — (T t {a,i(t)a sj(t')}) Their (Lehmann) spectral weight 

oo oo 

representations G^(E) = f du E +$-u an< ^ Ga(E) = I du E+i si p gj ^ E ^ s _ u yield either a symmetric imaginary part 

— oo — oo 

p(E) as spectral weight or the antisymmetric function pW(E) = sign(E) p(E) = — -Im{G(E)} with p( a >{E) = 
-p^(-E). 

We want to expand the antisymmetric function p^ a \E) in the complete orthonormal set of eigenfunctions of a ID 
harmonic oscillator. We shall describe below in detail how many states are necessary to describe with increasing 
accuracy the deviation of the numerically evaluated anti-symmetrized DOS p^ (E) from the (functional form of the) 
first excited state of the harmonic oscillator, ipi(x). We add one subsection in order to motivate the analogies between 
variables used in the DOS or P(h) and in the oscillator models' wavefunctions. 

Finally we shall compare the density of states of the fermionic SK-model with that of the Coulomb glass. Pioneering 
work by Davies, Lee, and Ricei^ and by Griinewald, Pohlmann, and Wiirtz 16 provided strong numerical evidence for 
the connection between the Efros-Shlovskii pseudogap and glassy order in Coulomb glasses, which intend to model 
the deep localized regime of disordered electrons interacting by Coulomb interaction. 

In these models, random chemical potentials were either modeled by a constant probability distribution on a finite 
width, while later on, Miiller and Pankov — used a Gaussian distribution to obtain an effective pseudo-spin model and 
analyzed the ordered phase in great detail^. 

Our objective is to display both the similarities and the differences between density of states of the fermionic SK- 
model and of the Coulomb glass. We report a new calculation on the replica symmetric approximation including the 
unstable temperature-regime and T = 0. We are then in a position to appreciate the different types of pseudogaps 
created by RSB. For 3D Miiller and Pankov found that the DOS vanishes quadratically at the pseudogap-center, while 
in the SK-case it behaves linearly for all dimensions. 



B. The oscillator models 



The one-dimensional harmonic oscillator is one of the most elementary models of quantum mechanics. Its eigenfunc- 
tions are usually written in terms of Hermite polynomials and Gaussian functions (while the 2nd linear independent 
solution of the differential equation, which is a confluent hypergeometric function, is excluded by boundary condi- 
tions). Using this complete orthonormal set of basis functions for representations of other models' solutions, does 
not imply a priori any physical relationship with the harmonic oscillator. Hence in general one cannot expect to gain 
much from such a representation. 

The apparent resemblance between P(h) of the T = SK-model (and hence of the density of states p{E) of its 
half-filled fermionic model extension^) on one hand, and the modulus of the first excited state wavefunction |?/;i(a;)| 
of a harmonic oscillator on the other hand, motivated us to look for the number of basis states needed to describe 
well the deviations. Another faint suspicion was feeded by the replica-symmetric spin glass solutions at T — 0, which 
has a split Gaussian shape and reminded (partially) of the oscillator ground state. 

Analyzing spin glass data in terms of such an expansion, we first found that less basis functions are needed than 
expected. Furthermore, for an excellent approximation (see below) it turned out that an unconventional shift was 
needed. This model can be described as 

H osc = ^ + -muj 2 (\x\ - x ) + c 3 \x\ 3 + c 5 \x\ 5 + ... (5) 

where p = —ih-^- The introduction of the nonanalytic shift helped to select only the first excited state in order to 
match accurately the DOS-data by only the first excited state. For simplicity we set again h = 1 and m = 1. 
The energy of the oscillator in the n-th excited state, given by 

E<£l = J dxTp*(x)H osc ip n (x), 

while the SK-energy E$k at T = will then be given by 

o n=0 

E SK =F SK (T = 0) = J dhhP(h) = J dE E p(E) (6) 

— oo — oo 

and using the results for P(h) or p(E) we confirm the recently derived high-precision value of the SK-energy at T = 
correctly up to 0(1O -8 ). 
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Anticipating the result let us represent the SK-energy E$k as an integral over the quantum oscillator wavefunction 

by 

o 

E SK = J dxx\Mx)\ (7) 

— oo 

which is the direct translation of Eq.©. 

The SK-energy Esk has little to do with the oscillator energy in the first excited state mentioned above. We do 
not introduce the special quantum oscillator as a simpler 'replacement '-model for the SK spin glass, but for the time 
being rather focus on the mere demonstration that its first excited state reproduces the functional form of P(h) or 
p(E), which are both two important quantities of the spin glass order at T = 0. 



C. Distinguishing differential equations of the T = SK-model from Burgers equation 

In order to motivate the corresponding sets of variables needed in the functional mapping of DOS, P(h) and one 
hand and the oscillator eigenfunctions on the other, we wish to return to the differential equations of the spin glasses 
at zero temperature. 

It was observed that the Parisi scheme led to a recursive relation for the so-called exponent-correction-function, 
denoted by expC, in arbitrary finite order of RSB (discrete Parisi scheme). The recursion relation turned into a 
partial differential equation 9,17 in the so-called continuum limit. 

Here, for the purpose of comparing it with a Burgers equation, we apply a different simple transformation. It is 
sufficient to start from the original equation as given in Refi^) for fields < h < oo (see also Ref iA£ which includes 
h = 0) 



—^-d a expC{a 7 h) + ~ 
q'(a) 2 



2a—d h j expC(a,h) 



a (dhexpC(a, h)Y 



= 0. 



(8) 



where the order function q(a) is known to be monotonic and its derivative q'(a) vanishes only in the limit a 
Applying the transformation 



Y(a, h) := expC(a, h) + \h\ 



da a q (a) 



(9) 



simplifies further the PDE. We had given a numerically satisfying analytic model function q(a) for T — in RefJ^, 
which differs everywhere within < a < oo by only less than O(10 -3 ) from the exact solution. This function is 
monotonic and we may switch to the unique inverse function a{q) and, even better, choose t = 1 — q as a new 
independent variable, in order to remove the inconvenient sign. We call this second variable pseudo-time, since, 
by means of the transformation ([5]), the PDE now matches a diffusion equation with diffusion constant | and one 
nonlinear term. Apart from the inverse order function a(r) it resembles a Burgers equation as discussed just below. 



dY(r,h) _ 1 
~~ 2 



-^Y(T,h) + a(r) 



dY(r, h) 
dh 



(10) 



Apparently we obtain a modified Burgers equation with a coefficient a(r) which diverges for small argument like 
the inverse square root for q — > 1, hence a(r) ~ r -1 / 2 . In recent papers we concluded that {a — 0,T = 0} and 
{a = oo,T = 0} are too different critical points of RSB hosted by the SK-model. The latter limit belongs to q — > 1, 
while the small-a limit corresponds to q ~ a — > and a(r = 1) = 0. 

In order to complete our discussion above, we recall that the Burgers equation d t u(t,x) = Dd%u(t,x) — 
u(t,x)d x u(t,x) turns into 

8 t <j>{t, x) = Dd 2 x <t>{t, x) - X {d x (/)(t, x) f 

by means of u — d x (j> and dtd x cj> = d x dt4>- 

It is well-known that the Burgers equation has a redundant nonlinearity, and that it can be transformed into a 
diffusion equation by means of the Cole-Hopf transformation^. 

In the differential equation for the SK-model, Eq. (|10[) . however, the coupling function a(r) of the nonlinear term 
prevents a linearization by means of the Cole-Hopf transformation. This may indicate that an exact analytical solution 
of the SK-model at T — is much harder than the exact KPZ-solution in one dimensio n 19 ' 20 . This argument does 
not yet refer to the particular initial condition of the SK-model at T — 0. 
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D. Randomly stirred SK-model differential equation versus KPZ-equation 

It is natural to consider also a (non-thermal) noise perturbation of the T = SK differential equation. In case of the 
pure Burgers equation this perturbation is called random stirring and results in the KPZ-equation, which describes 
random growth of a surface height above a substrate during particle deposition. 

Let us perturb the SK-energy by a non-thermal <5-correlated noise function ry(t, h) such that the characteristic 
differential equation ([TU]) takes a form comparable with the KPZ-equation 

du{t, h) ^ D ^ ^ x) + x (y ^ ( ^ ^2 + ^ x) 



dt 

dY(r,h) _ 1 
dr ~~ 2 



— Y(r,h)+a(r) 1 



dh 2 v ' ' w V dh 



T)(T,h) (12) 



The renormalization of the one-dimensional KPZ-equation yields the critical exponents exactly 2 ^ and it is known 
that simple rescaling did not lead to those exponents, except for the Edwards- Wilkinson behaviour near the trivial 
fixed point. This fixed point is however unstable for one space dimension. A simple power counting analysis which 
compared different terms of the differential equation had no access to the behavior near the stable finite fixed point 
for the effective coupling^. An exact solution of the 1+1 dimensional KPZ-equation with flat initial conditions was 
presented by Calabrese and Le Doussali^. 

Let us reconsider the question whether the modified PDE, which describes the SK-model with randomly perturbed 
energy, may however allow power counting to be successful. Its disadvantage, which is the divergent coefficient function 
a(r) ~ l/v^~ °f the nonlinear term, may render this possible. 



1. Scaling regimes in the randomly-stirred SK-model differential equation 

We have a good analytical model for the T = order function q(a) for all a and almost exact knowledge of 
the asymptotic behaviour close to the critical points at a — and a = oo, as derived from a scaling theoretical 
explanation of high-precision numerical data^. This can be used to estimate the role played by the different parts 
of the differential equation. The diffusive part dominates the small a-limit and the nonlinear term the large a limit, 
respectively. 

Under this assumption that the noise term does not essentially change the large-a property 1 — q ~ 1/a 2 and hence 
the coefficient function a(r), we can estimate by power counting the competition between terms of the differential 
equation and, in analogy with the treatment of the KPZ-equation 2 -, eventually obtain power laws and even critical 
exponents. 

Rescaling the pseudotime byr-> 6 z r, the analog of the real-space variable by h — > b h, and the field Y by Y — > b a Y 
one gets for the small r-regime (large a-regime or q(a) ~ 1 — .41/ a 2 in terms of the spin glass order function q(a)) 

dr 2W l ' )+ ^ \ dh J + V[ ' > 

In the large a-regime or small r-regime the nonlinear term dominates over the diffusive one and one gets 

3 ah 
a = — , z = 1 = — . 

4 2 2 

This simple power counting result does not replace a complete renormalization group calculation. Note that in this 
language of dimensionless pseudotimes r = 1 — q the above scaling behaviour refers to a critical short-time limit. It 
describes the behaviour of the function Y(q) or expC(q) for q- values close to the diagonal of the Parisi-matrix. 

For the long r-regime, which refers to small a in the original equation, one has q(a) ~ a behaviour instead and the 
diffusive term dominates over the nonlinear one. For this regime one is left with a randomly stirred diffusive equation 
in leading order and the behaviour is the equivalent of Edwards- Wilkinson (EW) behaviour in the KPZ growth-model 
for one space dimension (and one time dimension). In this case one gets simply z = 2 and a — \ as for the EW-limit 
of the KPZ-equation. 



III. EVALUATION OF THE ZERO TEMPERATURE DENSITY OF STATES 



We consider the spectral function of the fermionic SK model, which has the virtue of being mappable to the internal 
field distribution P(h) of the standard SK-modelii. Thus we can make use of the functional identity p(x) — P(x). 
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The density of states function p K (u) at a fixed order k of replica-symmetry breaking steps can be described by 
the following hierarchical integrals below. They form the basis for the present numerical evaluation of p K (E) for 
k = 0, 1, 2, 100 and on a dense grid of energies E. An RSB solution for small orders k < 4 can be found in22, which 
also considered arbitrary filling and an additional Hubbard interaction. 

The following formula for the DOS holds for arbitrary order k. Its evaluation allows to understand precisely the 
RSB-flow through the quasi-continuous regime at large k = O(10 2 ) towards k = oo. 



p K (E, X i) = - i exp(a 1 (\E\-xi))U\E\,0)6(\E\- Xl ) (13) 

27iy<?i - q_2 

which involves k nested integrals /„ (E, h) given by the recursion relation 

oo 

I n (E,h)= J dti gih-h'^q^-q^iD^ih'))^- 1 I n ^(E,h'). (14) 

— oo 

The initial condition is given by the free-propagator like forms 

E - 

~%i - 92) 
with 

9ix ' v) = ^^(-lj)- (16) 

A second recursive structure is furthermore needed in Eqs ll4l which now concerns the nested integrals of type D n (x). 
Their recursive structure is given by 



LAE.h) ■ .j-^, *' (15) 



exp( 



(h-xy 



D n {x)= I dh ^ (J n _ 1 (h))— (17) 

V5n+1 — q n +2 



with the initial condition 



D (x) =cxp[-aix+-a 1 (qi-q 2 ))-\l-erf(- 



2 7 2 \ y/2( qi - q 2 ) 7 

exp( ai x + \a\{ qi q 2 ))\ ( 1 - erf C X Z a(qi ~ q2) -)) ■ (18) 



2 iw 



V 2 (9i - 52) 



The parameter sets {a^}, {qi}, and xi: which are required in this recursive set of equations, depend on the order k. 
They are inferred from a previous calculation, which reported their values^ for all k up to the maximum order 200. 
They were obtained by cxtremization of the SK free-energy and their T = subset can be used since the fermionic 
SK- and the SK-model features coincide for T = and half- filling (p, = 0), as mentioned in the model section Hi Al In 
particular the non-equilibrium susceptibility \i vanishes for k — > oc~, which gives rise to the perfect pseudogap. 



IV. ANALYSIS OF THE NUMERICAL HIGH ORDER RSB DATA 



Numerical results for the T = density of states have been evaluated for all k = 1,2..., 100 orders of replica 
symmetry breaking. One observes an incredibly fast convergence such that after k = 10 any high k numeric solution 
can be effectively treated as the k = 00 result. For the following calculations the data for k = 100 are used. The quality 
has been checked by comparing pwo(E) with the fixed point function lim p K {E) evaluated by using Pade-series. The 

difference is negligibly small rendering either choice equally good. 

Finding an analytic form of the density of states function is not only interesting but also useful for further calculations 
(e.g. Greens-Function for the model). In the course of finding a good fit for the data, it turned out to be quite practical 
to antisymmetrize the data by hand, thus getting rid of the non-analytical point at the origin. 

Looking at this new data-set the similarity to the first excited state of the harmonic oscillator becomes rather 
striking. However the detailed analysis shows some necessary corrections as described in the next section. 
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p(E) P(E) 




FIG. 1: (left figure:) The plot shows the numerical result pioo(E) for the DOS at k = 100 (blue, symmetric curve) compared 

with its partial sum X^n=o c ioo 7 /'™(-£') °f the harmonic oscillator decomposition and the first excited state ipi(E) of the ap- 
proximate Hamiltonian H approx . In the lower right corner the difference Ap(E) = pioo(E) — p apP rox{E) is plotted for the two 
approximations. Note that both approximations are correct up to 0(1O~ 3 ) at each point. 

(right figure:) This example displays, for selected energies in the pseudogap regime, the flow of density of states data for the 
fermionic SK- Model under increased RSB-order k. The approach of the fixed point DOS p{E) = lim K ^ 00 p K (E) with slope 0.3" 
is obvious (even from a subset of the available k). 



A. Calculating the overlap of p(E) with harmonic oscillator eigenfunctions 

The eigenfunctions of the harmonic quantum oscillator, described by the Hamiltonian Ho = ^x 2 + in dimen- 
sionless variables, given by 

1 / a 

ip n {x) = H n (x)exp — 

i!2™ V 



with the Hermite polynomials H n (x), form the basis for our calculations. 
Calculating the overlap coefficients 

C W = J dx^ n {x)p^\x) 
quantifies the deviation of the function p( a ' (x) from \p\ (x) . 



B. Construction of an anharmonic Hamiltonian and functional mapping 

It turns out that it is possible to perturb the harmonic oscillator in a way that the DOS (here p^li 00 ) coincides with 
the functional form of the oscillator's first excited state. The perturbation used in our case is rather unusual, since 
it involves the absolute values of the £ operator, being \x\. Since the overlap-coefficients go to zero rather fast 
and the computation of the change in all eigenfunctions due to the perturbation is impossible, we restrict ourselves 
to a subset of the first 22 eigenfunctions. That renders an algebraic treatment of the Hamiltonian as a 22 x 22 matrix 
possible. 

We use the fact that eigenstates and eigenvalues of an anharmonic oscillator can be well approximated by diago- 
nalization of the non-diagonal matrix representation (ipk\H an h\ipi) m a subset of eigenstates \ip) of the unperturbed 
harmonic oscillator. 

H 0mn = (n + ^)S mn (19) 



dx ip m {x)\x\ip n {x) (20) 
tf p = ff + a|£|+/3|i| 3 +7|i| 5 (21) 
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In order to find the correct coefficients a, ft, 7 for the given perturbation it is useful to define the following function 
l(a, /3,7, A), introducing a new eigenvalue A. In order to find the correct variational parameters solving H p c = Ac it 
is useful to define the function 



Z(a, /3, 7, A) = \\H p [a, p,j]c- A g|| 



(22) 



which controls the numerical accuracy. The function is constructed in a way such that the correct coefficients minimize 
it. 

1. Anharmonic quantum oscillator eigenstate modeling the spectral function 

The results for the leading parameters a, /?, 7 of the oscillator potential and for the energy eigenvalue of the first 
excited state are given by 



a = -0.899165, /3 = -0.003268, 7 = 0.0000405, A = 0.364335 



(23) 



and the square distance from the minimum is I = 9.9 10 -7 . The minimum was found by employing the Minimize- 
function implemented in Mathematica. 

While the focus is on the first excited state, Fig. 3 includes the three lowest states for the given potential parameters 







\ 1.0 





0.5 



FIG. 2: The three lowest eigenstates of the nonanalytic cubic potential V(x) = ^x 2 + a\x\+ ft\x\ 3 +y\x\ 5 are displayed together. 
The position of their base-line indicates the corresponding energy eigenvalue. The first excited eigenstate (odd parity) reproduces 
the density of states p(x)sign(x) of the fermionic SK-model, and hence P(h)sign(h) of the standard SK-model too, with full 
Parisi symmetry breaking at zero temperature. 



Even the cubic term could be neglected and still the perturbed Hamiltonian produces a very good fit of the calculated 
DOS with slightly modified parameters 



a = -0.9225, = 7 = 0, A = 0.3444 



(24) 



and I = 7.3 10~ 6 . 

Thus, by keeping only linear and quadratic terms, one may rewrite the oscillator model with an unconventional 
quadratic potential shifted by a constant, but using \x\ instead of the normal space operator x. Alternatively, one 
may use the standard real space operator together with an unconventional shift, which flips its sign at x = 0. The 
potential function as included in Fig. 3 shows a double- well feature induced by the unusual shift; a difference between 
the cubic and the quadratic approximation is invisible on the given scale. 
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The resulting approximate Hamiltonian assumes the following form 



appiox — 



\ [f + (|z| - 



(25) 



In order to control the finite cutoff imposed on the number of harmonic oscillator basis functions, the effective model 
with nonanalytic shift and small anharmonic terms, as derived above, can be plugged into a Schrodinger equation, 
which is then solved numerically. 

Neglecting anharmonic terms of 0(|a:| 3 ) because of their smallness but keeping only the nonanalytic shift (hence 
the |a;|-term), the approximate model Hamiltonian -ff a pprox leads to the differential equation 



The numerical solution of this equation confirms the results of our large but finite matrix diagonalization approach. 
This means, not only the desired first excited state ipi is well approximated, but the other stationary states too. We 
used this as an independent control of our procedure. 

V. A COMPARISON WITH THE MEAN FIELD SOLUTION OF THE COULOMB GLASS 

Since we focus in the present article on characterization and classification of replica symmetry breaking at lowest 
temperatures including T = 0, we wish to compare its role in two different host models by means of pseudogap 
features. There seems to be enough evidence that the host-model, which undergoes this type of symmetry breaking, 
can have an influence. While the universal critical behavior of replica symmetry breaking at T = is primarily caused 
by its hierarchical structure, it can yet belong to different universality classes depending on the host-model. 

Recent work on the so-called Coulomb glass (CbGl) offers an interesting comparison with the fermionic SK-type of 
RSB. Both models can be viewed to describe the insulating phase of a disordered itinerant fermionic systems - the 
source of glassiness being either a random chemical potential (CbGl) or a frustrated spin interaction respectively. 

A selfconsistent structure was derived by Muller and Pankov for the Coulomb glass models. Their starting point, 
was the transformation to an effective spin glass model, justifying a mean field description by means of the long-range 
nature of the bare Coulomb interaction. The resulting model appeared as a type of SK-model with an effective 
1/r-dependent pseudo-spin interaction. This spatial dependence provokes a dimensional dependence. 

An intriguing feature of the Miiller-Pankov selfconsistent theory is the large depletion of the fermionic density 
of states already far above the glass transition, identified by the Almeida-Thouless temperature. The dip of p(E) 
apparently becomes very pronounced already above T g = Tat, and the authors have shown the changeover into the 
final Coulomb pseudogap by taking the full Parisi RSB into account. The similarity with the Efros-Shklovskii result 
was demonstrated. The evaluation close to T = indicated a quadratic decay of the DOS in 3D, in agreement with 
the ES-law. The authors also reported some results for the two-dimensional case. 

In order to see the quantitative genuine effect of replica symmetry breaking in the CbGl-model, we evaluated 
their self-consistent coupled equations without replica symmetry breaking also below T g . As initially suspected, 
the difference between solutions with broken and unbroken symmetry becomes larger with temperature decreasing 
towards zero. At first sight, the pseudogap seemed to be shaped already above T g (resembling a lot the curves of the 
fermionic SK-model in a certain range below its T c = T g , see Fig[3]), with only small corrections to be expected from 
symmetry breaking as T decreases further. But we found that, within the replica-symmetric self-consistent structure, 
the depletion starts to form a hardgap, which, under certain assumptions, can even be elaborated analytically. The 
numerical results down to T = 10 -3 and moderate disorder are shown in Fig|3l Using the same notation and units 
asi£, the hardgap half-width equals h , which is evaluated to be h = 0.894306684. This value could be given more 
accurately, but is is not needed here, it also shows that only parts of the calculation can be done numerically. 

Thus, one might say that the pseudogap-like form, in part developed already above T gi would be destroyed without 
RSB as T — > 0. Hence RSB is again responsible for the existence and shape of the pseudogap. This is one of the 
crucial points for a comparison with the fermionic SK-model: the replica symmetry breaking is fully responsible for 
the pseudogap in the low temperature limit. 

The second point is that the CbGl-pseudogap shows, according to Muller and PankovJ^, a d-dependent power law 
and in particular a quadratic decay near the center of the pseudogap for 3D, while the SK-model shows linear behavior 
(for all d). Thus the host of replica symmetry breaking can have an effect on its critical behavior, since, in our point 
of view, the small energy deviation from the gap-center belongs to the critical point a = oo of the order function q(a) 
as derived recently in Ref. 10 , and hence to small pseudo-times r = 1 — q. 



d%.vp(x) + (x 2 + 2a\x\ - 2X)ip(x) = 



(26) 



10 




-4 -2 ~ 2 4 -6-4-2 2 46 



FIG. 3: left panel shows density of states (DOS) of the fermionic SK-model for temperatures above T c (red) and below (blue), 
including the exact analytical result at T — (dashed-black) without replica symmetry breaking. Right panel: For comparison, 
the analogous replica-symmetric (RS) re-calculation on the basis of the Miiller-Pankov Coulomb glass (CbGl) equations is 
shown. The depletion of the DOS almost looks like a pseudogap already at temperatures far above the Almeida- Thouless 
T g but, as the temperature decreases further towards zero, a hard gap is formed in the RS-approximation of the CbGl as in 
the fermionic SK-model. Displayed curves are shown for temperatures above T g (red), at T g (green), below (blue) down to 
T = 10 -3 . The true pseudogaps (see above for the fermionic SK) and the Coulomb glass— obey different power laws and differ 
strongly from the shown RS-solutions. 

VI. CONCLUSIONS AND REMARKS 

The main purpose of the present article is to show that the internal magnetic field distribution and the density of 
states as central quantities of SK spin glass physics can be described by the functional form of the first excited state 
of a weakly anharmonic quantum oscillator with nonanalytic shift. This meant the harmonic potential to be of the 
form (x — Xq sign{x)) 2 = (\x\ — xq) 2 instead of x 2 . 

Puzzled by the unexpected relationship between spin glass and oscillator with a non- analytical shift, we searched 
for other singular anomalous oscillator models. We became aware of recent work by Ritort on wedge-shaped oscillator 
potentials^. He reported glassy behavior and aging effects in so-called generalized oscillator model(s) GOM, including 
the case of a pure and positive | x | -potential. 

In our independent study the |x|-term has opposite sign and required (at least) a; 2 -term for stability. The nonanalytic 
linear term was required as an optimized shift to reproduce correctly the p(E)-da,ta (and by analogy P{h) as well) 
with a minimal number of oscillator base functions. 

Let us also mention recent numerical work by Boettcher et al 24 on small- and moderate-sized spin glass models 
(including SK) for the magnetic field distribution at T = 0. 

The field distribution of the 3D Coulomb glass, and also those of XY- or Heisenberg spin glass models, show rather 
quadratic small field behavior. This raises the question whether they can be represented by oscillators too, but not 
by one single excited state. For continuous dimensions D one may also wonder how a p(E) ~ (E — Ep) ^ 1 can be 
obtained. Perhaps fractional derivatives can generalize the differential equations such that these power laws can be 
represented by oscillator models. These questions appear quite open and might lead to more insight into the spin 
glass oscillator relationship. 
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